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ABSTRACT 

We present the discovery of three z > 5 quasars in the AGN and Galaxy Evolution Survey (AGES) 
spectroscopic observations of the NOAO Deep Wide-Field Survey (NDWFS) Bootes Field. These 
quasars were selected as part of a larger Spitzer mid-infrared quasar sample with no selection based 
on optical colors. The highest redshift object, NDWFS J142516. 3+325409, z = 5.85, is the lowest- 
luminosity z > 5.8 quasar currently known. We compare mid-infrared techniques for identifying z > 5 
quasars to more traditional optical techniques and show that mid-infrared colors allow for selection 
of high-redshift quasars even at redshifts where quasars lie near the optical stellar locus and at z > 7 
where optical selection is impossible. Using the superb multi-wavelength coverage available in the 
NDWFS Bootes field, we construct the spectral energy distributions (SEDs) of high-redshift quasars 
from observed f?iy-band to 24/^m (rest-frame 600 A - 3.7/^m). We show that the three high-redshift 
quasars have quite similar SEDs, and the rest-frame composite SED of low-redshift quasars from the 
literature shows little evolution compared to our high-redshift objects. We compare the number of 
z > 5 quasars we have discovered to the expected number from published quasar luminosity functions. 
While analyses of the quasar luminosity function are tenuous based on only three objects, we find that 
a relatively steep luminosity function with \l/ tx L~ 32 provides the best agreement with the number 
of high-redshift quasars discovered in our survey. 

Subject headings: quasars: general; quasars: emission line 



1. INTRODUCTION 

Understanding the evolution of quasars through cosmic 
time allows us to study the history of accretion of matter 
onto supermassive black holes in the nuclei of galaxies. In 
addition, quasars at high redshift provide the lighthouses 
needed to probe the conditions under which galaxies ini- 
tially formed near the epoch of re-ionization. Samples 
of quasars at high redshift, however, are still quite small 
and most high-redshift quasar surveys are only sensitive 
to the most luminous examples of these objects. We 
know very little about the population of low-luminosity 
quasars at these early times. 

Surveys for low-redshift quasars have been quite 
successful and have led to quasar samples suit- 
able for detailed statistical studies of the evolu- 
tion of the quasar luminosity function to moder- 
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ate redshifts (ISchrriidt fc Oreenlll 9831 IBovle et a.lJ |200f : 
|Fa,n et a1Jl2001 rial IWolf et a,l . 1 120031 ICroorn et allB 
[ Richards et alJ 120051 iBarger et all 120051 iBrown et alJ 
120061 1 Jiang et alJl2006t iRichards et al.ll2006j) . Censuses 
of high-redshift objects have opened the Universe even 
back to the era of reionization, but the number of such 
objects is still small. Deep spectroscopic follow-up of i- 
dropout objects in the Sloan Digital Sky Survey (SDSS) 
( York et aTl 1200(1) have found 19 quasars at z > 5.7 



( Fan et al.ll200iai2lMl2lMl2liol . IZheng et all (200C 
and lChiu et alJ l)2005|) found 7 quasars at z ^ 5 based on 
SDSS photometry which comple ment the 17 quasars at 
5 < z < 5.4 found in SDSS DR3 iAbazaiian et alJl200l 
spectroscopy selected as outliers from the stellar locus 
(Schn eider et alJl2005() . SDSS imaging, however, is rel- 



atively shallow (zab < 20.5 mag) and thus the quasars 
found using SDSS photometry probe only the luminous 
tail of the quasar luminosity function. Furthermore, 
the use of optical photometry for selecting high-redshift 
quasars is biased against heavily reddened objects, as the 
observed optical flux from objects at z > 3 probes the 
rest-frame ultraviolet emission which is most heavily af- 
fected by dust and the strong absorption from neutral 
hydrogen in the intergalactic medium. 

The difficulty of conducting deep, wide-area multi- 
color surveys has meant that only a handf ul of lower- 
lumin osity sources have been reported. iStern et alJ 
(2000) found a faint quasar (M B = -22.7) at z = 5.50 
in a small area imaged approximately 4 magnitudes 
deeper than SDSS. A quasar at z = 5.189 (M B ~ 
—23.2) was selected based on 1 Ms of X-ray imagin g 
in the Chandra Deep F ield North ijBarger et al.1 12002^1 . 
iDiorgovski et alJ l|2003D imaged the field around the 
quasar SDSS J0338+0021 at z = 5.02 and discovered 
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a second quasar at z = 4.96, a magnitude f ainter than 
the SP SS detection limit (Mb = -25.2). iMahabal et all 
(2005) imaged the field around the z = 6.42 quasar 
SDSS J1148+5251 and report the discovery of a faint 
(M B = -24.3) z = 5.70 quasar. 

Unfortunately, the number of successful searches for 
faint high-redshift quasars is comparab le to the num- 
ber of surveys with negative results. iBarger et alJ l)2003f) 
found no additional high-redshift quasars in the full 2 
Ms imaging of the Chandra Deep Field North that were 
not de tected in 1 Ms of observation and iCristiani et al.l 
(2004) found no high-redshift objects based on deep X- 
ray imaging of the Hubb le Deep F i eld N orth and Chandra 
Deep Field South. iWillott et all (l2005l> imaged 3.8 deg 2 
of sky in i and z and lSharpTtTaTl p004|) imaged 1.8 deg 2 
in VIZ reaching 3 and 2 magnitudes deeper than SDSS, 
respectively, but neither survey detected any new low- 
luminosity, high-redshift quasars. 

In this paper, we present photometry and spectroscopy 
of three new high-redshift (z > 5) quasars using the 
multi-wavelength photometry ava ilable in the N OAO 
Deep Wide-Field Survey fNDWFS: IJannuzi fc DevtToM 
Dey et al., in prep.) Bootes field and spectroscopic ob- 
servations from the AGN and Galaxy Evolution Survey 
(AGES; Kochanek et al., in prep), which provides red- 
shifts for several highly-complete samples of quasars to 
Jab = 22 mag, nearly 2 mag fainter than SDSS selects 
high-redshift objects. With our three newly discovered 
z > 5 quasars, we compare the number density inferred 
from this survey to predictions based on quasar luminos- 
ity functions in the literature. 

This paper is organized as follows: in §2, we sum- 
marize all the multi-wavelength photometry used in 
the paper and in §3 we discuss our spectroscopic ob- 
servations. We compare the mid-infrared selection of 
quasars utilized in this work to optical criteria used in 
the past in §4. Finally, we place our survey in con- 
text with past studies of the quasar luminosity func- 
tion and consider future high-redshift quasar searches in 
§5. All optical photometry presented here are corrected 
for foreground galactic reddening using the dust maps 
of iSchlegel. Finkbeiner. fc Davis I dl9 98). We use AB 
magnitudes for all bands ijOke II1974|) . although the pho- 
tometric catalogs from the NDWFS and FLAMINGOS 
Extragalactic Survey (FLAMEX; lElston et all l2005|) 
present Vega magnitudes 1 . Flux measurements from 
Spitzer are converted to AB magnitudes using tuab = 
23.9 — 2.5 log(,/i,/l/iJy). Also, when quoting optical, 
near -infrared, and IRAC ph otometry, we use SExtrac- 
tor ijBertin fc Arnoutsl 1199^1 MAG.AUTO magn itudes 
(whic h are comparable to Kron total magnitudes; iKronl 
1980) due to their small systematic errors and uncertain- 
ties at faint fluxes. When calculating luminosities, we use 
a (O m ,f2A) = (0.3,0.7) flat cosmology and H = 70 km 
s -1 Mpc~ . 

2. MULTI- WAVELENGTH PHOTOMETRY 

2.1. NOAO Deep Wide- Field Survey 

We utilize the deep optical (BwRI) and near-infrared 
(K s ) imaging of the 9.3 deg 2 Bootes field provided by 
the third data release from the NOAO Deep Wide-Field 

Bw,AB = Bw,Vegsu, RAB = -Rvcga + 0.20, I AB = -?Vcga + 0.45, 

Jab = ^Voga + 0.91, and K StA B = Ks tVcga _ + 1.85 



Survey l|Jannuzi fc DevllT999|) . A full description of the 
observing strategy and data reduction will be presented 
elsewhere (Jannuzi et al., in prep; Dey et al., in prep) 
and the data can be obtained publicly from the NOAO 
Science Archive 2 . The NDWFS catalogs reach B\y,AB ~ 
26.5, Rab ~ 25.5, I A b ~ 25.3, and K S;AB ~ 23.2 at 50% 
completeness. 

2.2. zBootes 

We imaged 8.5 deg 2 of the sky, covering 7.7 deg 2 of the 
NDWFS Bootes field, between 29 Janua ry 2005 and 31 
March 2005 in the z'-band with 90Prime \ Williams et alJ 
I2004|) on the Bok 2.3m telescope on Kitt Peak. The 
90Prime imager offers a 1 deg 2 field of view when 
mounted at prime focus on the Bok telescope. Expo- 
sure times range from 1-2.5 hours throughout the field 
with typical seeing of l'/8. Images were flatficldcd using 
observations of the twilight sky. For each night of obser- 
vations, we stack all of the dithered z'-band images to 
produce a high signal-to-noise ratio image of the fringing 
pattern in the detector. This master fringe image is then 
scaled and subtracted from each object frame to remove 
the strong fringing pattern. 

We calibrate the astrometry and photometry of 
these z'-band images using public imaging from the 
Sloan Digital Sky Survey D R4 ijYork et al.l l200?i 
lAdelman-McCarthv et aTll2005[) . Cross comparisons be- 
tween zBootes and SDSS show a O'.'l rms dispersion in the 
astrometry and a 5% scatter in the photometry for bright 
stars (z < 19). This 5% scatter is likely a combination 
of photometric calibra tion errors in SPS S (expected to 
be on the order to 2%; llvezic et al.ll2004j) and imperfect 
fringe removal and flat fielding in the zBootes imaging. 
Furthermore, we find that objects observed in overlap- 
ping zBootes fields have an rms scatter of 5% in the final 
photometry. The astrometry between the NPWFS and 
zBootes are slightly offset (Aa « 0'.'4, A5 « O'.'l). When 
matching objects in each zBootes field to the NDWFS 
catalogs, we remove the local mean astromctric offsets in 
both right ascension and declination from zBootes posi- 
tions, resulting in astrometry that agrees to 0'.'2 rms. As 
the exposure times and observing conditions were vari- 
able throughout the survey field, the limiting depth of 
the catalog depends on location in the survey area. The 
typical 3cr depth is 22.5 mag for point sources in a 5 
arcsecond diameter aperture. Full details of the data re- 
duction and a full release of the z'-band imaging catalogs 
will be presented in a future data-release paper (Cool, in 
prep) . 

2.3. FLAMEX 

The FLAMINGO S Extragalactic Survey (FLAMEX) 
ijElston et alJ [2005) provided some of the J and K s 
photometry used in this work. The Florida Multi- 
object Imaging Near-IR Grism Observational Spectrom- 
eter (FLAMINGOS) on the Kitt Peak 2.1m telescope 
was used to image 4.7 deg 2 of the NDWFS Bootes field 
to a limiting depth of K St AB ~ 21.1 detecting approxi- 
mately 150,000 sources (5a). The FLAMEX catalogs are 
publicly available 3 . 

2 http://www.archive.noao.edu/ndwfs 
|http: / / www.noao.edu/noao/noaodeep | 

J http: //flamingos. astro. ufi.edu/cxtragalactic/overview. html 
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2.4. IRAC Shallow Survey 
The IRAC Shallow Survey <|Eisenhardt et al.l2004|) ob- 
served 8.5 deg 2 of the sky in the NDWFS Bootes re- 
gion a t 3.6, 4.5, 5.8, an d 8.0 /imwith the IRAC instru- 
ment ijFazio et all 12004) on Spitzer. This survey found 
« 270,000, 200,000, 27,000, and 26,000 sources brighter 
than 5cr limits of 12.3, 15.4, 76, and 76 /iJy (correspond- 
ing to limits of 21.2, 20.9, 19.2 and 19.2 AB mag) in 
each of the four IRAC bands l|Eisenhardt et alJ 120041 
IStern et al.ll2005l . 

2.5. MIPS Imaging 

The NDWFS Bootes field was also observed at 24, 
70, and 160 /im wi th the Multiband Imaging Pho- 
tometer for Spitzer l|Rieke et alJ 120041) as part of the 
Spitzer IRS team's Guar anteed Time Observing pro- 
grams l|Houck et al.ll200"5|) ; only the 24 /im photometry 
is considered here. The 24 /im imaging covers 8.22 deg 
of the NDWFS Bootes field and reaches a 80% complete- 
ness limit of 0.3 mJy l)Brown et alJ l2006). 

3. SPECTROSCOPIC OBSERVATIONS 
3.1. AGN and Galaxy Evolution Survey 

The AGN and Galaxy Evolution Survey (AGES, 
Kochanek et al., in prep) has obtained complete spec- 
troscopic samples of galaxies and quasars using several 
multi-wavelength selection techniques in the NDWFS 
Bootes field. Spectra of « 20,000 objects were taken with 
Hectospec, a 300 fiber robotic spectrog r aph on the MMT 
6.5m telescope l)Fabricant et all 119981 120051: iRoll et al.l 
1998). Data reduction was completed using HSRED, 
a modified version of the SDSS spectroscopic pipeline. 
Dome flat spectra were used to correct for the high- 
frequency flat-field variations and fringing in the CCD, 
and, when available, twilight sky spectra provided a low- 
frequency flat-field correction for each fiber. Each Hec- 
tospec configuration has approximately 30 fibers dedi- 
cated to measuring the sky spectrum which is subtracted 
from each object spectrum. Simultaneous observations 
of F-type stars in each configura tion are cross correlated 
against a grid of Kurucz models l)Kuruczl|1993j) to derive 
a sensitivity function for each observation, thus linking 
the observed counts to absolute flux units. The final re- 
duced spectra cover the 3700A to 9200A spectral range 
at resolution R « 1000. 

AGES quasar target selection occurred in two stages. 
In the first phase of the project, AGES obtained red- 
shifts of nearly all point sources with Rab < 21.7 mag 
that were either X-ray sou rces with 4 or mor e counts 
in the the XBootes survey ijMurrav et al.ll2005l) or were 
MIPS 24/im sources brighter than lmJy with non-stellar 
24/tmto J-band colors. This led to a sample of roughly 
900 spectroscopically identified AGNs as well as red- 
shifts of nearly 9000 galaxies selected using a variety 
of other techniques. By combining the AG ES galaxy 
and AGN redshift samples. IStern et alJ l)2005f) confirmed 
the ease with which quasars could be recognized using 
mid-infrared colors. In the second phase of AGES spec- 
troscopy, we extended the optical flux limit to Iab = 22 
and added IRAC color selection to our methods of iden- 
tifying quasars. In detail, mid-infrared selected objects 
considered here were required to have [3.6] — [4.5] > —0.1 
and either [3.6] < 20.8 or a detection at 5.8 or 8.0/imand 



be classified as point-sources in the NDWFS imaging. In 
AGES, targets were classified as point sources if their 
SExtractor CLASS_STAR parameter was larger than 0.8 
for at least one of the Bw, R, or /-bands. The AGES 
mid-infrared selection cut used for point sources dis- 
penses with the [5.8] — [8.0] color restriction given in 
IStern et al.l 1)200 5?) which is intended to minimize the con- 
tribution of z > 1 galaxies which are not common to our 
optical flux limits. While we only consider the sample of 
point-sources selected in the mid-infrared from AGES in 
this paper, it is worth noting that when selecting AGN 
from extended source in AGES, the vertical selection cuts 
from IStern et alJ l|2005l) are imposed to reduce contami- 
nation of star-forming galaxies on the AGN sample. The 
shift from an i?-band to an /-band spectroscopic limit led 
to a significant increase in the number of z > 4 quasars 
in the AGES database, including three z > 5 objects. 
The second phase of AGES target selection also included 
a 24/im -selected quasar sample. These objects were re- 
quired to be point sources, have /24/im > 0.3mJy, and 
have /-band 3'.'0 diameter aperture magnitudes, I^,ab, 
such that I^ab > 18.45 — 2.51og(/ 2 4 Mm /mJy). This I- 
band selection criterion rejects normal stars from the 
sample but retains AGNs. Two of the three z > 5 quasars 
selected based upon their IRAC colors were also selected 
based upon their 24/im photometry. 

3.2. Auxiliary Observations 

In order to augment the AGES discovery spectra, spec- 
troscopic observations of two of the high-redshift ob- 
jects discovered in this work were completed with the 
Mufti- Aperture Red Spectrograph (MARS) on the 4-m 
telescope on Kitt Peak and the DEIMOS spectrograph 
l)Faber et a l. 2003) on the Keck-II telescope under poor 
observing conditions. These data were reduced using 
standard techniques. Also, to achieve higher quality 
data, two of the objects described in this paper were ob- 
served multiple times within AGES itself. Coaddition of 
all of the available spectra for each object was performed 
to improve the final signal-to-noise ratio of the data. 

4. RESULTS 

We have identified three new high-redshift (z > 5) 
quasars observed as part of the Iab < 22 AGES spec- 
troscopic sample. Figures and El show the /-band 
finding chart and spectra for each of the three high- 
redshift objects. Tables 1 and 2 list the observed pho- 
tometry and rcdshifts of the three new high-redshift 
quasars and Table 3 shows the derived rest-frame B- 
band luminosity and the luminosity of each object at 
1450 A. In calculating the rest-frame luminosities, we 
assume a power-law SED of the form f v cx v a " where 
a u = —0.5. The spectra shown in Figure |3 are the 
coadded spectra from all the observations of each object: 
NDWFS J142516.3+325409 was observed three times in 
AGES as well as with MARS and DEIMOS, NDWFS 
J142937.9+330416 was observed once in AGES and with 
MARS, and NDWFS J142729.7+352209 was observed 
twice within AGES. Each of the discovered objects show 
the clear signature of a broad asymmetric Lya emission 
line. As the signal-to-noise ratios of the final spectra 
are only modest, redshifts for each object were deter- 
mined based upon the break in the Lya line profile rather 
than on weak emission lines which are poorly detected. 
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TABLE 1 

Optical and Near-Infrared Photometry 



Object Name 


Redshift 


B w 


R 


I 


z' 


J 


K s 


NDWFS J142937.9+330416 
NDWFS J142729.7+352209 
NDWFS J142516.3+325409 


5.39 
5.53 
5.85 


> 27.1 

> 26.3 

> 25.8 


22.84 ±0.09 
24.20 ±0.21 
23.99 ± 0.11 


21.24 ±0.05 
21.85 ±0.07 
21.57 ±0.06 


21.02 ±0.07 
21.93 ±0.11 
20.68 ±0.06 


20.40 ±0.12 


20.81 ±0.19 
20.41 ±0.19 



All magnitudes are on the AB system 



TABLE 2 
Spitzer Photometry 



Object Name 


Redshift 


[3.6 fim] 


[4.5 fim] 


[5. 




[8.0 Aim] 


[24 fim] 


NDWFS J142937.9+330416 


5.39 


20.81 ±0.08 


20.24 ±0.09 


> 


19.2 


19.65 ±0.13 


17.31 ±0.15 


NDWFS J142729.7+352209 


5.53 


20.27 ±0.06 


19.84 ±0.07 


> 


19.2 


19.88 ±0.14 


17.53 ±0.18 


NDWFS J142516.3+325409 


5.85 


20.36 ±0.06 


19.91 ±0.07 


> 


19.2 


20.74 ±0.18 


> 17.7 



All magnitudes are on the AB system 



NDWFS J142937. 9 + 330416 NDWFS Jl 42729.7 + 352209 



< -20 
-40 
-60 
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FlG. 1. — Finding charts for the three z > 5 quasars discovered 
in the AGES survey. The images are drawn from the NDWFS /-band 
imaging of the Bootes region and are each 2 arcminutes on a side. 
North is up and East is left. The location of the quasar is marked by 
vertical lines in each panel. 



Quoted redshifts thus have errors of order Az ~ 0.02. 

The quasars included in this paper were not selected 
with standard dropout techniques which identify ob- 
jects with very red optical colors indicative of a strong 
spectral break. These techniques have been success- 
ful in identifyi n g quasars at redsh i fts a b ove 5 (e.g. 
Fan et all l2000T l2001d 120031 12001 l2005t IStern et all 
20001 iDiorgovski et alJl2003t iMahabal et alJl2005|) . but 
have several weaknesses. At z ^ 5.5, quasars cross the 
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FlG. 2. — Final stacked spectra for each of the high-rcdshift quasars 
discovered by our work. These spectra represent the coadded spectra 
for the repeated observations with Hcctospcc as well as observations 
with MARS and DEIMOS. Two prominent quasar emission lines, Lya 
and NV, are marked with vertical dashed lines and the measured red- 
shift of each object is listed in the upper left corner of each panel. 
All three of the quasars show the strong Lya emission line truncated 
by neutral hydrogen absorption blueward of the emission line center 
characteristic of high-rcdshift quasars. 



stellar locus in optical color-color space making opti- 
cal selection problematic. The inclusion of near-infrared 
photometry can break this degeneracy to some extent, 
but obtaining this data can be time consuming. 

The high-redshift quasars studied here were all se- 
lected on the basis of their mid-inf rared colors using a 
scheme similar to that presented in IStern et alJ |2005). 
In brief, quasars, with their roughly power-law spectra, 
have redder mid-infrared colors than z < 1 galaxies. The 
IRAC color-color space for AGES objects shown in Fig- 
ure |3| illustrates the separation of the high-redshift ob- 
jects discovered here from stars and low-redshift galax- 
ies. We also indicate the criteria used for AGN selec- 
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tion given bv iStern et al.1 ([2005) (dashed line) and the 
AGES point-source mid-infrared selection criterion (dot- 
dashed line). Extended source s are target e d as AGN 
in AGES if they meet the full IStern et all lEffll cri- 
teria. Low-redshift [z < 0.5) objects, which are pre- 
dominately extended galaxies, are shown with x-marks 
while objects with z > 1 are marked with circles and are 
dominated by unresolved AGN. At [5.8] - [8.0] ^ 1, the 
locus of low-redshift galaxies clearly crosses the AGES 
point-sourc e mid-infrar e d sele ction criterion but remains 
outside the IStern et alJ l)2005f) selection region, illustrat- 
ing the need for the separate point-source and extended- 
source selection criteria used by AGES. The three high- 
rcdshift quasars studied here (diamonds with error-bars) 
are well-separated from low-redshift galaxies throughout 
this color-color space. 
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FlG. 3. — The mid-infrared colors of objects detected in the AGES 
spectroscopic survey. The symbol type and color denote the rcdshift 
of each object. The x-marks show the 0.01 < z < 0.5 objects in 
AGES. This sample is primarily composed of galaxies, although some 
low-redshift AGNs arc also included. The empty circles (red) show 
objects with z > 1.0 and arc dominated by AGNs. The triangles (blue) 
mark the location of stars in this color-color space. The filled diamonds 
(green) show the location of the z > 5 qua sars in this color-co lor space. 
The dashed line shows the region used by Stern et al. I 2005) to select 
AGNs. The vertical lines are used to protect against the presence 
of z > 1 galaxies in the AGN sample, but as galaxies at z > 1 are 
fainter than our optical flux limits, these cuts are not used in the AGES 
mid-infrared point-source quasar selection shown by the dot-dashed 
line. When selecting AGN from extended objects in AGES, the vertical 
lines are used to reduce contamination from low-redshift star-forming 
galaxies. 

One obvious question arises due to our mid-infrared 
selection of high-redshift quasars: does the population of 
mid-infrared selected quasars differ from quasars selected 
using more traditional optical techniques? The spectra 
in Figure [3 though of only modest signal-to- noise ratio, 
show no obvious differences from objects presented in the 
literature and would suggest no strong difference between 
high-redshift quasars selected based upon their optical or 
mid-infrared colors. Figure [1| shows the I— z versus R — I 
color-color space for point sources observed in AGES as 



well as th e color track of high - redsh ift quasars predicted 
using the IVanden Berk et a.lJ ll200~ft SDSS quasar com- 
posite modified by the additio n of absorption blueward o f 
Lya using the prescription of iSongaila fc Cowid |2002). 
The three quasars discovered here reside in the region 
traditionally used to select quasars with z > 4.5 (e.g. 
iRichards et al.ll2002|) further indicating that the quasars 
selected in this study are similar to those selected in 
the past. Both NDWFS J142516. 3+325409 and NDWFS 
J142937. 9+330416 lie quite close to the stellar locus in 
these colors, however, illustrating one of the main diffi- 
culties of traditional optical color selection techniques. 
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FlG. 4. — Optical color-color space for AGES point sources. The 
symbols arc defined as in Figure^] Additionally, the dashed line shows 
the expected colo rs of quasars at z > 4.8 predicted using the SDSS 
quasar composite IVanden Berk et all2001D modifie d by the addition of 
absorp tion blueward of Lya using the prescription of Soneaila & c Cowiel 
12002). The three high-redshift quasars discovered in this work occupy 
the same location in color-color space used in the past to optically- 
select z > 4.5 quasars. At the same time, two of the three objects lie 
very near the stellar locus in this color-color space, illustrating one of 
the main difficulties in using optical colors alone to select high-redshift 
quasars. 



Figure shows the spectral energy distribution from 
the observed Bvv-band to 24 /im (approximately rest- 
frame 600 A - 3.7 ^m) for each of the high-redshift 
quasars identified in this work. Again, all three of 
these objects were selected based upon their mid- 
infrared colors; the two lowest-redshift quasars were 
also selected based on their 24/im fluxes. For compar- 
ison, we have also plotted the average quasar SEDs 
from two studies in the literat ure; the top line shows 
theE vis. Risaliti. fc Zamoranil l)2002|) radio-quiet quasar 
composite while the lower line shows th e average SEP 
of ffpite er-observed SDSS quasars from IRichards et alJ 
(2006b]). The broad photometric properties of the high- 
redshift quasars are quite similar to each other as well 
as to the low-redshift composite spectra. The main dif- 
ference between the low- and high-redshift objects is 
the lack of absorption due to neutral hydrogen at the 
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TABLE 3 

Luminosities of Discovered Quasars 



Object Name 


Redshift 


M B 


A/1450 


NDWFS J142937.9+330416 


5.39 


-26.00 


-25.52 


NDWFS J142729. 7+352209 


5.53 


-25.15 


-24.67 


NDWFS J142516. 3+325409 


5.85 


-26.52 


-26.03 



highest frequencies in the composite spectra. The high- 
est redshift object, NDWFS J142516.3+325409, shows a 
strong break between observed 4.5/miand 8.0/im, but 
this is likely due to poor signal-to-noise ratio in the 
8.0/im photometry as this object is near the flux limit of 
the 8.0/im imaging. The rest-frame UV to optical color of 
NDWFS J142729.7+352209 is redder than the other two 
quasars studied here as well as both of the comparison 
composites, possibly indicating enhanced dust extinction 
in this object. Compared to a sample of 58 quasars at 
3 < z < 3.5 from AGES, the rest-frame UV to optical 
color of NDWFS J142729.7+352209 is not unique; it falls 
well within the distribution of colors of lower-redshift ob- 
jects. 
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FlG. 5. — Broad-band spectral energy distributions of the three 
high-redshift quasars presented here. The quasars are ordered by 
redshift with the highest-rcdshift object on the top. The top two 
lines show two composite quasar spectral energy distributions for 
comparison; the top line shows the radio-quie t quasar SED from 
lElvis. Risaliti. fe Zamoranil I2002D while the bottom comparison shows 
the com posite of Spitzer observed SDSS quasars from iRichards et al.l 
I 2006b). The high-redshift SEDs arc very similar to each other and the 
low-rcdshift composites arc broadly similar to the SEDs of the high- 
redshift objects although the SEDs of the three z > 5 o bjects presented 
here have much stronger breaks in the blue from absorption from neu- 
tral hydrogen. The large ratio between the 4.5/imand 8.0/imflux of 
the highest-redshift object, NDWFS J142516. 3+325409 is likely a re- 
sult of low signal-to-noise ratio photometry as this object is close to 
the detection limit in that band. 

Figure shows the optical versus mid-infrared color- 
color space for AGES point sources. In this color space, 



all three of the discovered z > 5 quasars are separated 
from the locus of low-redshift galaxies, Galactic stars, 
and z < 5 quasars. We show one possible selection 
method for high-redshift quasars by the dot-dashed line 
in FigureEl Within this region, the AGES source catalog 
contains 19 point sources with I < 22 mag; 14 of these 
were spectroscopically observed. Of these 14 targets, 9 
objects are stars, 2 are lower-redshift (z = 2.11,3.57) 
quasars, and 3 have z > 5. Assuming the 4 objects with- 
out redshifts are not located at high redshift, this selec- 
tion results in a minimum efficiency of ~ 17%, though 
with only 3 quasars, this efficiency measurement is rather 
uncertain. The large contamination from stars arises due 
to increasing photometric errors in the [3.6] — [4.5] colors 
as the stars approach the IRAC flux limit and thus much 
of the contamination could be mitigated with a deeper 
IRAC imaging survey. 

5. DISCUSSION 

The quasars reported here add three new low- 
luminosity, high-redshift quasars to the slowly growing 
catalog of these objects. All three of these objects are 
fainter than quasars at similar redshifts found in the 
SDSS, and NDWFS J142516. 3+325409 is the lowest lu- 
minosity z > 5.8 quasar currently known. The number 
of these objects should grow quickly as the next gener- 
ation of deep wide-area surveys are completed, opening 
the door to understanding the nature of low-luminosity 
quasars near the epoch of reionization. 

As the number of low-luminosity z > 5 quasars is still 
small, the details of the luminosity function o f these ob- 
j ects is poorly constrained. Extrapolating the iFan et aT] 
Il2001bft quasar luminosity function (QLF), determined 
using 3.5 < z < 5 quasars with —27.5 < M1450 < —25.5, 
to higher redshifts and lower luminosities, we would ex- 
pect a density of 0.28 deg -2 to the optical flux limit of 
the AGES spectroscopy or 2.2 quasars with z > 5 in 
the AGES survey area for a complete, optically limited, 
survey. 

In order to estimate the effect of the IRAC flux lim- 
its on the expected number of quasars in our survey, 
we create a sample of 58 quasars at 3 < z < 3.5 from 
the full AGES catalog with existing FLAMEX photom- 
etry. These objects include quasars selected using the 
full suite of AGES selection techniques including X-ray, 
radio, or 24/im fluxes or their optical or mid-infrared col- 
ors and thus likely exhibit a broad range of broad-band 
photometric properties. We augment this sample with 
17 quasars in the sam e redshift range listed in the SDSS 
Third Quasar catalog ijSchneider et al.ll2005j) which were 
also detected in 2MASS. For each object in this sample, 
we convert the observed R — K s (or r — K s for SDSS 
quasars) to a rest-frame ultraviolet-to-optical spectral 
slope, ajjOj by assuming the measured broad-band col- 
ors are the result of a pure power-law SED with /„ oc v a . 
Next, we create a m o ck cata log of high-redshift quasars 
using the IFan et all (|2001b|) QLF to assign each mock 
object a redshift and UV luminosity at 1450 A. The ob- 
served distribution of UV-to-optical spectral slopes de- 
fines the distribution of fc-corrections that are applied to 
the mock catalog to convert the rest-frame luminosity at 
1450A to the observed flux at 3.6/xmfor each object. We 
then use the observed distribution of mid-in frared spec- 
tral indices, cumir — 0.73±0. 84. measured bv lStern et all 
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(2005), to assign each mock quasar a flux at 4.5, 5.8, and 
8.0 /Ltm . Finally, we apply the AGES mid-infrared flux 
limits and selection criteria to measure the fraction of 
high-redshift objects which are missed by the AGES mid- 
infrared selection criteria. We find that approximately 
35% of the mock quasars that pass the I < 22 mag op- 
tical flux limits are missed when the IRAC flux limits 
arc included. It may be a concern that the addition of 
the SDSS quasars, which, unlike the AGES sample, were 
all selected based on their optical colors, may bias this 
estimate; we have verified that omitting these objects 
from the calculation does not affect the final complete- 
ness more than a few percent. It should be noted that 
the above calculation assumes quasars at z > 5 have the 
same broad-band photometric properties as z ~ 3 ob- 
jects; this assumption can only be tested after a large 
sample of high-redshift quasars with multi-wavelength 
photometry is collected. In the I < 22 AGES sample 
of 2153 mid-infrared selected objects, 83% of the targets 
obtained valid redshifts. Combining this 83% spectro- 
scopic completeness with the effects of the mid-infrared 
flux limits, we estimate that our overall completeness is 
54% and thus we would expect to find only 1.2 quasars 
at z > 5 for a QLF with a slope of W oc L -2 6 . To then 
find three quasars is somewhat unlikely, since the Pois- 
son probability of finding 3 or more objects when 1.2 are 
expected is only 12%. 
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FlG. 6. — Optical versus mid- infrared color-color space of AGES 
spcctroscopically observed point sources. The symbols and colors are 
as in Figure The dashed line illustrates the color track of quasars 
above z ~ 4.8. The z > 5 quasars are well-separated from the loci of 
low-rcdshift galaxies, stars, and lower-redshift quasars. The region de- 
fined by the dot-dashed line illustrates one way to select high-redshift 
quasars by combining optical and mid-infrared photometry. The min- 
imum efficiency for selecting z > 5 quasars in this region is 17% with 
the largest portion of the cont amination arising due to stars with low 
signal-to-noisc IRAC photometry which scatter to redder [3.6] — [4.5] 
colors. Deeper mid-infrared imaging would increase the color measure- 
ment accuracy, significantly reducing the stellar contamination and in- 
creasing the quasar selection efficiency. 

Based on a sample of 12 quasars at z > 5.8. lFan et alJ 



(2004) found that the high-redshift QLF was best fit 
by a steeper QLF (\& oc L -3 ' 2 ) than that measured by 
IFan et alJ t200lU) . If we instead use the best fit b right- 
end slope and normalization from IFan et al.l l)2004|) but 
keep the redshi ft evolution of the quasar number density 
determined bv IFan et alJ l|2001b|) . we would expect to 
find 6.2 quasars at z > 5 in a complete optically-limited 
survey of the AGES area. This value drops to 3.5 when 
the effects of our mid-infrared flux limits and spectro- 
scopic incompleteness are added. This value agrees quite 
well with the 3 objects found in our survey, but, at this 
point, this agreement is only suggestive due to the small 
number of quasars in our sampl e and the s mall sample 
size and luminosity rang e used bvlFan et al.l ll2004|) to de- 
termine this QLF slope. lRichardsetHi!Ti(2006|) found the 
QLF slope to evolve from j3 = —3.1 (where "J oc L' 3 ) at 
z < 2.4 to progressively flatter slopes at higher redshift 
such that at z ~ 5, one would predict <; — 2.4. If the 
QLF at 5 < z < 6 is, indeed , as st eep at L~ 3 , the trend 
measured bv lRichards et afl 1)20061) must break down be- 
yond z ~ 4.5, but a larger sample of low-luminosity ob- 
jects is required to place any robust constraints on the 
evolution of the QLF to z ~ 5. 

We have shown that mid-infrared selection of quasars 
can be useful in the search for high-redshift objects to 
z ~ 6, thus avoiding some of the key problems inherent in 
optically based high-redshift quasar searches. At z <; 7, 
optical selection of quasars is impossible as the Lya emis- 
sion from objects at these redshifts is shifted out of the 
z'-band. In the near-infrared, the existing large area sur- 
veys (primarily 2MASS) are too shallow to detect even 
the brightest high-redshift quasars found by the SDSS. In 
the near future, the UKIDSS survey will image 4000 deg 2 
of the sky to 3 mag deeper than 2MA SS and should find 
abou t 10 quasars with 5.8 < z < 7.2 ijWarren fc Hewettl 
l20f)l . While near-infrared spectroscopy is ultimately 
the only means by which one can ensure any candidate 
is located at high redshift, the use of mid-infrared colors 
in the selection of these objects may provide a valuable 
tool in separating stars from quasars at the highest red- 
shifts. Figured illustrates the expected [3.6] — [4.5] color 
of quasars to z ~ 8 based upon the locally determined 
SPSS quasar composite spectrum ijVanden Berk et all 
I2001[K We also show the colors predicted from a power- 
law [a v — —0.5) continuum with Balmer emission lines 
having the same strength and width of those measured 
on the SDSS composite spectrum. 

As the strong emission from the Ha emission line red- 
shifts into the 4.5/imband, the colors of quasars red- 
den quickly near z ~ 5. Correspondingly, the sudden 
drop in color at z ~ 6.5 occurs as the strong Ha emis- 
sion line redshifts out of the 4.5/imband; by z ~ 7, the 
colors of high-redshift quasars redden as H/3 redshifts 
from the 3.6//mband into the 4.5/imIRAC bandpass. 
The overall trend toward bluer colors exhibited by the 
SDSS composite across this redshift range occurs as the 
IRAC passbands probe the transition between a spectral 
slope of a„ = —0.5 measured blueward of 5000 A and 
of a u = —2.5 redward of this wavelength in the SDSS 
composite spectrum. This transition likely occurs due 
to contamination of the AGN light by the host-galaxy 
in the low-redshift AGNs used to construct the com- 
posite spectrum. This contamination makes the colors 
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FlG. 7. — [3.6] — [4.5] color of quasars t o very high redshlft based 
on the SDSS quasar composite (solid line, |Vandcn Berk ct ail 120011) 
or a a — —0.5 power-law combined with the Balmcr emission lines 
observed in the SDSS composite spectrum (dash-dot-dot-dot line). The 
color cut imposed by the mid-infrared selection of quasars in AGES is 
shown by the dot-dashed line. Longward of rest-frame 5000A, the SDSS 
composite is dominated by low-luminosity AGNs and is likely heavily 
contaminated by host-galaxy light, making the colors slightly redder. 
For comparison, the location of the three z > 5 quasars discovered 
here are marked and the dotted lines on the right mark the colors of 
pure power-law spectra with different values of a u . The rapid rise at 
z — 5 occurs as the Ha emission line passes from the 3.6//mband into 
the 4. 5/j-m bandpass. Similarly, as the Ha emission line redshifts out 
of the 4. 5pm bandpass near z ~ 6.5, the colors drop quickly; the colors 
become redder again near z ~ 7 as the H/3 enters the 4. 5pm bandpass. 
Even at very high redshifts, the colors of quasars predicted by both 
tracks fall within the region used to mid-infrared select quasars in our 
survey. 



predicted using the iVanden Berk et all l)200H) compos- 
ite generally redder than those expected using the sim- 
ple model consisting only of a power-law continuum and 
Balmer emission lines. For both tracks, the mid-infrared 
colors of very high-rcdshift quasars remain redder than 



stars ([3.6] — [4.5] ~ —0.5) and low-redshift galaxies 
throughout the redshift range shown. As the large num- 
ber of deep, wide-area, Spitzer fields with corresponding 
deep optical and near-infrared photometry become avail- 
able, the ability to build large samples of high-redshift 
quasars, including low-luminosity objects, will allow us, 
for the first time, to probe the statistics of low- luminosity 
quasars at high redshift. 
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